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Granular particles vibrated in a fluid have been found to exhibit self-organization with attractive 
and repulsive interactions between the particles. These interactions have been attributed to the 
steady streaming flow around oscillating particles. Here we examine the steady streaming flow 
surrounding a vertically oscillating sphere using three dimensional particle tracking velocimetry. We 
present measurements of the flow with the sphere far from boundaries, close to the bottom wall of the 
tank, and near another oscillating sphere. The steady velocity field is found to disagree with available 
analytic calculations. When the sphere is oscillated near the bottom wall the entire topology of 
the flow changes, resulting in a larger repulsive region than expected. Previous experiments saw 
attraction between particles in the region where the flow around a single particle is repulsive. We 
conclude that advection in the streaming flow due to a single particle cannot explain the observed 
attractive and repulsive interactions, rather non-linear interactions between the flows around two 
or more spheres must be responsible. 



I. INTRODUCTION 

Steady streaming flows due to the oscillatory motion of 
a fluid are responsible for a wide range of fluid phenom- 
ena. Work on this problem goes back to Rayleigh^. Ap- 
plications of steady streaming flows include driving fluid 
with ultrasonic beams^, microfluidic transport!^, vesicle 
license"^, and streaming flows in the ear^. 

Recent experiments have observed self-organization of 
granular particles when they are vibrated verticall}*^ 
or horizontally^ in a fluid. Streaming flows were iden- 
tifled as the mechanism for self-organization, but there 
remained puzzling questions about the structure of the 
streaming flows and the attractive and repulsive interac- 
tions between particles. In this paper we present exper- 
imental 3D particle tracking measurements of the flow 
around vertically oscillating spheres. We show that for 
a streaming flow far from boundaries the available ana- 
lytic solutions do not agree with our measurements, likely 
because the theories are for limiting cases that are far 
from the experimental parameters. Measurements of the 
streaming flow at different distances from a solid bot- 
tom wall reveal that the wall changes the topology of the 
streaming flow. We also present measurements of the 
flow around two oscillating spheres. 

II. PARAMETER SPACE FOR OSCILLATING 
SPHERES 

In Figure [1] we show the parameter space relevant for 
a single oscillating sphere far from any boundaries. The 
vertical axis is the oscillation amplitude (half peak to 
peak) nondimensionalized by the particle radius, 

e = A/a . (1) 

The horizontal axis is the Reynolds number, 



Analytical solutions of the steady streaming flow around 
an oscillating sphere have three limits, which are indi- 
cated by Roman numerals in Figure [H The theoretical 
calculations are usually described in terms of alternate 
nondimensional paramters 
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Rcs = e Re = — — . 

The quantity Sqsc = \/^V^ thickness of the os- 

cillatory boundary layer. In Fig. [H solid lines are drawn 
at = 1 (upward sloping) and Res = 1 (downward 
sloping). The analytical solutions in limits I and II were 
found by Riley^^. Both these solutions are in the limit 
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FIG. 1: Parameter space for a single oscillating sphere in an 
infinite fluid. Plus symbols indicate data from this study; 
asterisk indicate data from Voth et. al^; triangles indicate 
data from Blackburn^. The Roman numerals indicate the 
limits described in Section |ll] for which calculations of the 
steady stream flow field are available. 
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of small Res (so that the streaming flow is dominated 
by viscosity) and small e (so that the oscillation ampli- 
tude is inflnitesimal) . The solution in limit I requires 
<C; 1 so that the oscillatory boundary layer thickness 
is much larger than the size of the sphere. In this limit 
the streaming flow is inward along the oscillation axis 
and outward at the equator. The solution in limit II re- 
quires 1 so that the oscillatory boundary layer is 
much smaller than the sphere size. Here the streaming 
flow has an inner recirculation zone and an oppositely di- 
rected outer recirculation zone that flows inward toward 
the equator and outward along the oscillation axis (see 
Van Dyke^^ or Fig. [3] for images of this flow). Solution III 
was obtained by Brenner and Stone^, who assume that 
the outer streaming flow is potential flow {Reg ^ 1 and 
e <^ 1). The flow fleld here is qualitatively similar to 
Solution II. 

The symbols in Fig. [1] show the parameters of the ex- 
periments in this paper, Voth et. al^, and the simulations 
by Blackburn^. Despite agreement between Solution III 
and experiments on spheres bouncing above a vibrating 
plate^, we will see that the existing theories are inade- 
quate to quantitatively describe the streaming flows in 
our experiments. Factors that contribute to the discrep- 
ancy include the non-inflnitesimal oscillation amplitude, 
the intermediate Reynolds number of the streaming flow, 
the effects of the bottom wall, and non-linear interactions 
between the flow around two spheres. 



III. EXPERIMENTAL SETUP 

A diagram of the experiment is shown in Figure [2l The 
sphere is held in the center of an octagonal plexiglas tank. 
The tank has an inner diameter of 30 cm and a height 
of 20 cm. The octagonal shape was chosen because it 
nearly has cylindrical symmetry, but still provides eight 
flat walls for optical access. The sphere is attached to a 
thin rod that moves on a linear motion rail. The vertical 
oscillations are driven by a stepper motor as shown in 
Fig. [Sfb). We primarily used a 1.91 cm diameter sphere 
attached to a 0.16 cm diameter rod. For some experi- 
ments we used a larger 5.09 cm sphere on a 0.33 cm rod. 
The vertical position of the sphere and the amplitude of 
the stroke could be adjusted continuously. The tank was 
fllled with a mixture of glycerol and water, adjusted to 
give the required viscosity. 

Three-dimensional particle tracking (3DPTV) was 
used to reconstruct the three dimensional positions and 
the velocities of the particles in the flow from stereoscopic 
camera images^^. 3DPTV has become a standard tech- 
nique in many areas of fluid measurement, particularly 
in the study of turbulent flows^^^^^^^^. As we show here, 
it can also be used to rapidly map the full three dimen- 
sional flow fleld of simpler flows. Compared with PIV, 
3DPTV requires lower seeding density, so it has lower 
spatial resolution. For periodic flows like ours, phase av- 
eraging can be used to obtain high spatial resolution for 
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FIG. 2: (a) A top view of the experimental setup drawn to 
scale. Two cameras at 90° from each other provide a stereo- 
scopic view of a volume of the the flow around the oscillating 
sphere. The overlap of the cameras is indicated by the quadri- 
lateral. Three theater lights provide symmetric illumination 
for both cameras, (b) A side view of the setup not drawn to 
scale. The circular motion of the disk sitting on a stepper- 
motor is converted to a linear up-and-down motion by the use 
of a connecting rod and a linear motion rail. The distance be- 
tween the lowermost point of the sphere during its stroke and 
the bottom plate will be referred to as d. 



the full three dimensional flow fleld from a single data set. 
We have implemented a simple 3DPTV setup using two 
1.3 Megapixel high-speed cameras for stereoscopic imag- 
ing. This limits the seeding density to signiflcantly less 
than is possible with 3 or 4 cameras^^, but the setup was 
suflicient for our needs. Each camera has 4 GB of RAM 
which allowed 3200 images to be acquired continuously. 
The cameras imaged a volume surrounding the sphere 
roughly 8 x 8 x 8cm^ as shown in Fig. [21(a). The tracer 
particles were alumino-silicate microspheres 80-120 jam in 
diameter with density 0.7-0.9 g/cm^ (supplied by Trelle- 
borg Fillite, Inc.). Three 750 Watt theater lights gave 
symmetric illumination. 
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To find the three dimensional coordinates of the tracer 
particles, we match the image of a particle on one camera 
to the image of the same particle on the other camera. 
Each position on a camera image corresponds to a parti- 
cle that exists somewhere along a ray in space. The three 
dimensional position of the particle is determined by find- 
ing rays from the two cameras that intersect. Converting 
image positions to rays in space requires a camera model 
and an accurate calibration to determine the parameters 
of the model. We used a simple distortion free camera 
model with seven parameters defining camera location, 
viewing direction, rotation and magnification. To ob- 
tain the calibration parameters, we imaged a calibration 
mask in the glycerol mixture and measured the real space 
positions of each dot. A non-linear fit of these measured 
positions to the calibration model gives the necessary cal- 
ibration parameters for that camera. The system could 
determine the three dimensional position of a 100 /im 
particle with an accuracy of about 10 jam. 

Since we are primarily interested in the streaming flow 
created by the oscillatory motion of spheres, we typi- 
cally take one image per period of the sphere oscillation. 
Tracking particles through these phase locked images al- 
lows direct measurement of the steady flow. The cameras 
were triggered by a photogate aligned with 16 uniformly 
spaced holes around the circular disk. Usually we cov- 
ered all but one hole to resolve only the streaming flow. 
For some experiments we took 16 images per period to 
observe the oscillatory flow, i.e. how the recirculation 
zones changed during one driving cycle. 

Experimental difficulties limited our access to the full 
parameter space of our problem. Firstly, the rod which 
is attached to and drives the sphere had some influence 
on the flow fleld. It produced an up/down asymme- 
try in the streaming flow even when the sphere was far 
from the wall, where up /down symmetry is usually ex- 
pected. This effect became less important as Rcg was 
increased. Secondly, the continuous input of heat from 
the theater lights caused some convection in the ffuid. We 
installed infrared mirrors and absorbing glass, but where 
the streaming ffow is very weak (at small Re^ small e, or 
far from the sphere) the turbulent convective ffow is still 
visible in the data. Finally, at high glycerol concentra- 
tions (for many runs we used 92% glycerol) we found that 
the index of refraction of the ffuid became large enough 
that polystyrene particles, which we had planned to use 
as tracer particles, became invisible. The replacement 
particles had a density of {0.7 — 0.9 g/cm^) which was sig- 
niffcantly less than that of the ffuid (1.24^/cm^), so the 
particles would rise very slowly due to buoyancy. This 
had a negligible effect on the measured velocities, but 
after a long experimental run, there would be no tracer 
particles left in the region around the sphere. All three 
of these problems were less signiffcant for strong steam- 
ing ffows at higher Re and e. This is why our data is at 
higher Reynolds numbers than the data of Voth et. al^ 
and why we did not measure smaller e to more directly 
compare with the analytical results (see Fig. [1]). 
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FIG. 3: Streaming flow pattern in the center of the tank at 
Re = 30 and A/ a = 1.4. The plot shows height vs. radial 
position and the color indicates the radial velocity component, 
Vr. The position and size of the sphere are indicated by the 
semicircle. 



IV. EXPERIMENTAL RESULTS 

We ffrst looked at the ffow ffeld around a single sphere 
oscillating in the center of the tank, where no bound- 
aries should affect the ffow. We went on to study the 
changes that result from the presence of the bottom plate, 
where we observed a dramatic change in the ffow topol- 
ogy. The last set of experiments was conducted with a 
pair of spheres oscillating in phase. 



A. Flow field around a single oscillating sphere 

An example of the streaming ffow around an oscillat- 
ing sphere in the center of the tank is shown in Fig. [3l 
Images are taken once per oscillation at the center of the 
stroke, and then three dimensional positions of the tracer 
particles are determined. The cylindrical symmetry of 
the problem allows us to present the entire ffow ffeld by 
combining particles at all azimuthal angles, and plot- 
ting their height, vs. the distance from the oscillation 
axis, r. Plots of v^^ and horizontal cuts conffrm that the 
motion is almost entirely in the r — z plane in our ex- 
periments. (There are weak variations in the azimuthal 
direction, a result of convection currents, which lead to 
crossing tracks in regions of slow streaming ffow in the 
z-r-plots.) The radial velocity component Vr is shown in 
color. A few streamlines have been drawn in to indicate 
the direction of ffow. The particle tracking was optimized 
to observe the outer recirculations which cause the at- 
tractive interactions between particles. These ffows are 
much weaker than the inner recirculations, so this data 
does not resolve the inner recirculation zones. There are 
some mismatched particles resulting in a background of 
short spurious tracks, for example those inside the sphere. 
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FIG. 4: Steady inflow velocity Vrvs. r in a small stripe around 
z = 0. The scattered points are experimental data (Re = 30, 
e = 1.4, Rs = 42, |M|^ = 21.4), and are compared to the 
theoretical calculations of Brenner/Stone (dashed line) and 
Riley (sohd hne; NOTE: divided by 50 to fit on the plot). 
Background points are bad stereoscopic matches and should 
be ignored. 



1. Streaming velocity in the equatorial plane 

Figure [H shows the measured radial streaming veloc- 
ity as a function of radius for particles within 0.5cm of 
the equatorial plane. The flow is attractive at large dis- 
tances with a stagnation point 2cm from the center of 
the sphere. This figure also shows two theoretical calcu- 
lations of the streaming velocity. Riley^^ predicts that 
flows in limit II (from Fig. [1]) have 
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whereas Brenner and Stone^ predict that flows in limit 
III have 
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The agreement of the rate of attraction of two spheres 
with Eq. [5]was a result from Voth et. al^. 

Figure [H reveals a strong discrepancy between the the- 
ories and our experiment. The Brenner and Stone results 
(dashed line) has the right order of magnitude, but a sig- 
nificantly different functional form. Riley's result (solid 
line) has been divided by 50 to make it fit on the graph. 
The reason for the discrepancy is partly the impossibility 
of realizing an infinitesimal amplitude in the experiment. 
We think that the Brenner and Stone result does better 
because the experiment at Rcs = 42 is closer to the high 
Reg limit than it is to the infinitessimal Res limit used 
by Riley. However, viscous effects are likely still impor- 
tant in this Rcs range. Clearly, the inflow velocity in 
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FIG. 5: Steady inflow velocity Vr as a function of Re at fixed 
e = 0.7. A comparison between the experiment (filled shapes) 
and the Brenner/Stone-theory (open shapes) measured at 4a 
(circles), 6a (triangles), and 8a (squares) from the center of 
the sphere. 



the range r = a — 5a is not adequately captured by the 
available theories, and this is the range relevant for the 
experiments in Voth et. al^. Blackburn^ has simulated 
the flow around oscillating spheres and presents similar 
profiles of cycle-average radial velocities. His results also 
disagree with both of the theoretical predictions, but ap- 
pear to be consistent with our experiments. However, a 
direct comparison is not possible as his choices of oscil- 
lation amplitude, e = 0.1 and 10, are inaccessible to the 
experiment. 

The dependence of the equatorial velocity on Re is 
shown in Figure [5l We plot Vr at 4a (four sphere radii 
away from the center of the sphere), 6a, and 8a, and 
use e = 0.7 (half that of Fig. [4j). At each distance, 
the inflow velocity increases with increasing Re. This 
is in disagreement with Eq. [5] which predicts the veloc- 
ity should decrease with increasing Re as Ur/ (Aco) = 
— 0.53(i?e/2e)~-^/^(a/r)^. This discrepancy may be a re- 
sult of viscous effects on the streaming flow which are ig- 
nored in the potential flow solution used to obtain Eq. [5l 
Another feature of the inflow velocity is that it does in- 
crease with increasing oscillation amplitude e at fixed 
Rei^, but does not closely match either the linear de- 
pendence of Eq. [4]or the square root dependence above. 



2. Position of the equatorial stagnation point 

The stagnation point, where the inner and outer re- 
circulation zones meet in the equatorial plane, can be 
seen in Figures [3] and [H The position of the stagnation 
point, Tg, as a function of the Reynolds number Re is 
shown in Figure [6l Several runs are shown with different 
amplitudes and different sphere diameters. As Re is in- 
creasing, the stagnation point moves closer to the sphere, 
which is expected since the osci llato ry boundary layer 
thickness is decreasing: 5osc = \AV^ — <^(2e)"'^/^i?e~"'^/^. 
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ner which counteracts the increasing thickness of the os- 
cihatory boundary layer. An analytic solution for the 
streaming flow with non-infinitesimal oscillation ampli- 
tude would be very helpful in further understanding of 
these issues. 



B. Changes of the flow fleld due to the bottom 
plate 



FIG. 6: Position of the stagnation point as a function of Re 
at different A/ a. Filled diamond: A/a = 0.16, 2a = 2cm; filled 
triangle: A/a = 0.3, 2a = 5cm; open triangle: A/a = 0.5, 
2a = 5cm; filled circle: A/a = 0.7, 2a = 2cm; open diamond: 
A/ a — 1.4, 2a — 2cm; open square: A/ a — 2, 2a — 2cm. 
Note that rs is measured from the center of the sphere, so the 
surface of the sphere is at 0.5 on the graph. 
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The next set of experiments brought the sphere closer 
to the bottom plate of the tank. We took 5 single se- 
quences each having a different distance, from the bot- 
tom plate (see Fig. [2j) with fixed Reynolds number and 
amplitude {Re = 30, e = 0.7). Figure [8] shows that not 
only was the fiow geometry distorted by the presence of 
the boundary, but the topology also completely changed. 
The lower of the two far field recirculation zones suc- 
cessively closes in, gets compressed and eventually joins 
the upper of the two inner recirculation zones. Earlier we 
have used the position of the stagnation point to separate 
the inner repulsive recirculation and the outer attractive 
fiow. Here, this stagnation point has disappeared. An 
animation of the phase dependence of the fiow topology 
for Fig. [Hl^d) is available onlineii. 



FIG. 7: Position of the stagnation point rs as a function 
of A/ a at fixed Re = 25. The triangles are one oscilla- 
tory boundary layer thickness, Sosc, beyond the surface of 
the sphere. The surface of the sphere is located at 0.5. 



The position of the stagnation point already reveals a 
problem with previous interpretations of the attractive 
interactions between spheres. In Voth et. al^ the spheres 
attract to contact at approximately Re < 7 and e < 1.4. 
Our data, at Re = 20 and e = 1.4, shows that the stagna- 
tion point is 1.7 diameters from the center of the sphere. 
At lower Re, the stagnation point should extend farther. 
We conclude that the repulsive fiow inside the stagnation 
point extends across the distances where attraction was 
observed in Ref.^. Therefore, something more than sim- 
ple advection by the single particle streaming fiow must 
be involved in the fiuid mediated attraction. 

The dependence of the position of the stagnation point 
on oscillation amplitude is shown in Fig. [3 The size of 
the two inner recirculation zones has a maximum just 
below A/ a = 1, which is in contrast to the monotonic 
increase in the oscillatory boundary thickness. The rea- 
sons for this non-monotonic behavior are not fully under- 
stood. As the steady streaming velocity increases with 
increasing e, the steady boundary layer may become thin- 



C. Implications for interpreting fluid mediated 
interactions between oscillating particles 

The streaming fiow in Fig. [8] gives a vantage point 
from which we can evaluate the role of streaming fiows in 
the fiuid mediated particle interactions reported in Refer- 
ences^^^. Because of the topology change that unites one 
of the outer recirculations with an inner recirculation, 
the repulsive part of the the fiow extends much farther 
from the sphere than it does when the wall is not present. 
Even without the wall, the repulsive fiow extended too 
far from the sphere to explain the attraction observed in 
Voth et. al^. Now near the bottom wall, the repulsive re- 
gion extends to r = 4a and covers most of the attractive 
range. The current data is at higher Reynolds number 
{Re = 30 compared with Re 10), but Fig. [6] shows that 
the extent of the repulsive recirculation extends even far- 
ther at lower Re. The failure of Eq. [5] to describe the 
attractive fiow velocity further indicates that the attrac- 
tion and repulsion between vertically vibrated particles 
cannot be explained purely by advecting the particles in 
the single particle fiow fields. We conclude that the inter- 
actions previously observed^'— are the result of nonlinear 
interactions between the fiows around multiple particles. 
This conclusion has also been drawn from the study of 
the horizontally vibrated system by Klotsa et. al— . 
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FIG. 8: Changes to the streaming flow geometry as the sphere 
oscihates closer to the bottom wall of the container with fixed 
e = 0.7 and Re — 30. The plot shows height vs. radial 
position and the color indicates the radial velocity component, 
Vr. For (a) d/2a = 2.0, (b)= 1.5, (c)= 0.5, and (d)= 0.0. For 
each image the center of the sphere is at z = 0, and the sphere 
is at the center of the rod stroke. In (d) the sphere is almost 
touching the bottom wall at the bottom of its stroke. The 
position of the bottom wall is indicated by the black line. 
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FIG. 9: A periodic displacement plot in a 1-cm thick layer 
with Re = 30 and A/a = 0.7. This x = 0-layer is -0.5 <x < 
0.5. The plot shows y vs z and the color indicates the height 
velocity component, Vz. The position and size of the spheres 
are indicated by the semicircles. 



D. Flow field around 2 spheres 

Since non-linear interactions between the flows around 
two spheres must be responsible for the attractive and 
repulsive interactions, we have measured the flow around 
two oscillating spheres. The flow is now fully three di- 
mensional, which makes full use of the capabilities of 
3DPTV. We mounted two spheres on separate vertical 
rods and vibrated them in phase with a fixed distance of 
roughly 1 sphere diameter, 2a = 1.91cm, between them. 
Figure [9] shows the flow around two oscillating spheres 
near the bottom of the tank. Since this flow is not az- 
imuthally symmetric we plot only a slice of the velocity 
field that contains the vertical plane going through the 
centers of both spheres. Producing this image required 
10 times more data than the single particle flow fields 
since it could not be azimuthally averaged. The flow ge- 
ometry on the left and right edges of Fig. [9] shows roughly 
the same topology as the case of one sphere near the bot- 
tom wall (see Fig. [8]). Observation of the full flow field 
reveals that particles in the equatorial recirculation zone 
experience azimuthal drift toward the center where they 
are ejected in the upward jet. A 3D visualization of some 
of these trajectories is available onlinell. 



V. CONCLUSIONS 

We have demonstrated that a simple implementation 
of three dimensional particle tracking can provide accu- 
rate measurements of the flow fields around oscillating 
spheres. We set out to deepen our understanding of the 
fluid mediated attraction and repulsion between particles 
oscillating in a fluid, building on work by Voth et al^. 
This earlier work had observed clustering of particles os- 
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ciliated in a fluid, and found that for some parameters 
there could be repulsion and intriguing dynamic states. 
It also proposed that the attraction could be described 
by an analytical calculation of the streaming flow around 
single spheres. 

For parameters typical of the experimental observa- 
tions of fluid mediated interactions, our measured ve- 
locity flelds are quite different from available analytical 
calculations. This discrepancy is probably due to the 
theories being derived for limiting cases that do not ade- 
quately match experimental parameters. When a sphere 
is oscillating perpendicular to a flat plate, the stream- 
ing flow undergoes a dramatic change in topology due to 
the coalescence of two recirculation zones. This expands 
the range of the repulsive streaming flow near the sphere. 
Together these conclusions demonstrate that previous ex- 
planations^ gave oversimplifled descriptions of the flows 
that lead to the interactions between spheres. Some parts 



of the interactions can be qualitatively understood by 
single particle streaming flows. However, two particle ef- 
fects are important in controlling the interactions of par- 
ticles within a few diameters of each other at the stud- 
ied parameters. We have presented measurements of the 
streaming flow around a pair of oscillating spheres, but 
more work remains to be done in understanding how the 
non-linear interaction of the oscillatory flows leads to ob- 
served interactions between oscillating spheres. 
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